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Abstract. Effects of partonic interactions in heavy ion collisions at RHIC are 
studied in a multiphase transport model (AMPT) that includes both initial 
partonic and final hadronic interactions. It is found that a large parton scat- 
tering cross section is needed to understand the measured elliptic flow of pions 
and two-pion correlation function. 
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1. Introduction 

In heavy ion collisions at RHIC, a quark-gluon plasma is expected to be created 
in the initial stage. Many observables have been suggested as possible signals for 
the quark-gluon plasma. These include enhanced production of dileptons [ [l| and 
strange hadrons [ ^ , suppression of J /ip production [ ^) , Hanbury-Brown-Twiss in- 
terferometry between pions [ 0), elliptic flow of various particles [ ||, quench of jet 
production [ ^| , and fluctuations in the net charge [ |tJ , net baryon [ || as well as 
total baryon numbers [ ||. However, to establish these observables as unambiguous 
signals for the quark-gluon plasma requires that they are not affected appreciably 
by the final hadronic dynamics in the collisions. To study quantitatively these sig- 



nals, we have developed a multiphase transport model [ 10 that includes both the 
dynamics of initial partonic and final hadronic matters. With this model, we have 
studied not only the rapidity and transverse momentum distributions of various 
charged particles [ |ll| but also their elliptic flows [ [l2|]. Furthermore, we have in- 
vestigated multistrange baryon production [ J /if) suppression [ [b§ and pion 
interferometry [ [l5| in these collisions. We find that the model describes very well 
the experimental data on the charged particle rapidity distributions and transverse 
momentum spectra. However, these results are more sensitive to the hadronic in- 
teractions than the partonic interactions. On the other hand, both the charged 
particle elliptic flow and the two-pion correlation function are affected by partonic 
interactions, and to explain the observed large elliptic flow and the measured pion 
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correlation function in the experiments requires a large parton-parton scattering 
cross section in the initial partonic matter. 

2. Multiphase Transport Model 

In the AMPT model, the initial conditions are obtained from the HIJING model [ 
|l6| by using a Woods-Saxon radial shape for the colliding nuclei and including the 
nuclear shadowing effect on parton production via the gluon recombination mech- 



anism of Mueller-Qiu [ 17 . After the colliding nuclei pass through each other, the 
Gyulassy-Wang model [ |18fl is then used to generate the initial space-time informa- 
tion of partons. In the default HIJING model, these minijet partons are allowed 
to lose energy via the gluon splitting mechanism and transfer their energies to 
the nearby strings associated with initial soft interactions. Such a jet quenching 
is replaced in the AMPT model by explicitly taking into account parton-parton 
collisions via Zhang's Parton Cascade (ZPC) [ [lij. At present, only gluon elastic 
scatterings with a default cross section of 3 mb are included, so the partons do 
not suffer any inelastic energy loss as they traverse the dense matter. After par- 
tons stop interacting, they combine with their parent strings and, after an average 
proper formation time of 0.7 fm/c, are converted to hadrons using the Lund string 
fragmentation model [ ^| . Dynamics of the resulting hadronic matter is then de- 
scribed by a relativistic transport model (ART) [ ||lj, which has been improved 
to include baryon-antibaryon production from meson-meson interactions and their 
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The parameters in the AMPT model are determined by fitting the experimental 
data from Pb+Pb collisions at center-of-mass energy y/s — 17 'A GeV at SPS [ 
p3| . Specifically, to describe the measured net baryon rapidity distribution, we 
have included in the Lund string fragmentation model the popcorn mechanism for 
baryon-antibaryon production with equal probabilities for baryon-meson-antibaryon 
and baryon-antibaryon configurations. Also, to account for both the pion yield and 
the enhanced kaon yield in this reaction, we have adjusted the two parameters in 
the momentum splitting function used in the Lund string fragmentation model, 
following the expectation that the string tension is modified in the dense matter 
formed in the initial stage of heavy ion collisions. 

Since the probability for minijet production is very small in heavy ion collisions 
at SPS energies, the partonic stage does not play any role in these collisions. We 
find that final-state hadronic scatterings reduce the proton and antiproton yields, 
but increase kaon and antikaon production by about 20%. In contrast, the kaon 
yield in the default HIJING model is smaller than our final result by about 40%. 

3. Rapidity and Transverse Momentum Distributions 

Results on the charged particle rapidity distributions in central Au+Au collisions 
at RHIC are shown in the left four panels of Fig. [j] for center-of-mass energies of 
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56 A GeV (dashed curves) and 130 A GeV (solid curves) together with data from the 
PHOBOS collaboration [ ^4). The measured total charged particle multiplicities 
at mid-pseudorapidity at both energies are well reproduced by our model. Also, 
the p/p ratio predicted by the AMPT model is consistent with that measured in 
experiments [ ^5|. The transverse momentum spectra of pions, kaons, and protons 
in Au+Au collisions at ^/s — 130A GeV are shown in the right panel of Fig. [I]. 
The results from the the AMPT model are seen to also agree with the experimental 
data from the PHENIX collaboration [|6). 
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Fig. 1. Left four panels: Rapidity distributions of total and negatively charged 
particles, protons and antiprotons, charged pions, and charged kaons in Au+Au 
collisions at ^fs = h&A and 130A GeV. Right panel: Transverse momentum spectra 
of pions, kaons, and protons in Au+Au collisions at yfs = 130A GeV. The experi- 
mental data are for 5% most central collisions from the PHOBOS collaboration [ |24| 
for rapidity distributions and from the PHENIX collaboration [ |26| for transverse 
momentum spectra. Results from the AMPT model shown by curves are for impact 
parameters of b < 3 fm. 



Furthermore, the AMPT model describes very well the centrality dependence of 
the charged particle rapidity distribution measured by the BRAHMS collaboration 
[ as shown in Fig. [2] by dashed curves. Compared to the default HIJING model 
given by solid curves, hadronic scatterings are responsible for the broader rapidity 
distributions seen in the experiments. 

Without hadronic interactions, there is a significant increase in the numbers 
of total charged particles, pions, protons, and antiprotons at midrapidity. The 
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Fig. 2. Charged particle rapidity distributions for centralities of 0-5%, 5-10%, 20- 
30%, and 40-50% in Au+Au collisions at \fs — 130A GeV. The experimental data 
are from the BRAHMS collaboration [ |j . Theoretical results from the AMTP 
model are shown by dashed curves while those from the default HIJING model are 
given by solid curves. The figure is taken from Ref. [ I27J . 



kaon number is, on the other hand, reduced slightly. As a result, the ratios of p/p 
and K + /tt + in the absence of final-state hadronic interactions are 0.80 and 0.13, 
respectively, instead of 0.66 and 0.18 from the default AMPT model, which are 
more consistent with the data. Although the default HIJING [ ^8| gives a total 
charged particle multiplicity at midrapidity that is consistent with the PHOBOS 
data, including hadronic scatterings would reduce its prediction appreciably. 

Effects of partonic dynamics on final hadronic observables can be studied in 
the AMPT model by turning off the partonic cascade. We find that this leads to 
less than ~ 5% change in the final charged particle yields in Au+Au collisions at 
i/s = 130A GeV. The multiplicity distribution of hadrons is thus not very sensitive 
to parton elastic scatterings. 

4. Elliptic Flow 

Elliptic flow in heavy ion collisions measures the anisotropy of particle momentum 
distributions in the plane perpendicular to the beam direction. It results from 
the initial spatial anisotropy in non-central collisions and is thus sensitive to the 
properties of the hot dense matter formed during the initial stage of heavy ion 
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collisions [ [|, pSJ . In transport models based only on the parton cascade, the 
elliptic flow has been shown to be sensitive to the parton scattering cross section, 
and a large value can be obtained with a large cross section [ . 
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Fig. 3. Charged particle multiplicity (left panel) and transverse momentum (right 
panel) dependence of elliptic flow in Au+Au collisions at ^/s = 130A GeV. Data 
from the STAR collaboration [ |l| are shown by filled circles, while theoretical 
results for different partonic dynamics are given by curves. 

Results on the charged particle elliptic flow from the AMPT model for Au+Au 
collisions at center-of-mass energy of 130 A GeV are shown in Fig. |as a function of 
both the charged particle multiplicity (left panel) and transverse momentum (right 
panel). Theoretical values from the default AMPT model (no string melting) arc 
much smaller than that measured in the experiments by the STAR collaboration [ 
pl| . This is mainly due to the small fraction of energy that is carried by minijet 
partons and the lack of transverse collective motion of the strings. As a result, the 
elliptic flow is significantly reduced after minijet partons combine with strings and 
fragment to hadrons. 

Since the initial energy density of produced matter in central Au+Au collisions 
at RHIC is more than an order of magnitude higher than the critical energy density 
(~ 1 GeV/fm 3 ) expected for the formation of the quark-gluon plasma, keeping 
strings in the high energy density region thus underestimates the partonic effects 
in these collisions. To model the melting of strings into partons in regions of high 
energy density, we have extended the AMPT model to allow the initially excited 
strings to fragment into hadrons using the LUND fragmentation model and then 
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convert these hadrons to their constituent quarks but with bare masses, i.e., a 
meson is converted to a quark and an anti-quark, while a baryon is converted to 
three quarks. We further assume that quarks are produced isotropically in the rest 
frame of a hadron and start to interact only after a proper formation time given by 
the inverse of the hadron transverse mass. Scatterings among these quarks are then 
treated using the parton cascade ZPC. After the quarks stop interacting, we model 
the hadronization by combining the nearest quark and antiquark into a meson and 
three quarks into a baryon with the same flavor. The resulting hadrons are given 
an additional formation time of 0.7 fm/c in their rest frame and then imported to 
the ART hadronic transport model to take into account their rescatterings. 

With strings converted to partons, the initial energy originally stored in strings 
also contributes to the parton dynamics. This leads to a larger elliptic flow and also 
a stronger dependence on the parton scattering cross section, making it possible 
to determine the strength of partonic interactions from the final elliptic flow. As 
shown in Fig. ||, a large parton scattering cross section of 6-10 mb is needed to 
account for the observed dependence of charged particle elliptic flow on the total 
charged particle multiplicity and their transverse momentum. 

5. Pion Interferometry 

Particle interferometry based on the Hanbury-Brown Twiss (HBT) effect can pro- 
vide information not only on the spatial extent of an emission source but also on its 
expansion velocity and emission duration [}2| Ej3|. In particular, the long emis- 
sion time as a result of the quark-gluon plasma to hadronic matter phase transition 
in the initial stage of heavy ion collisions is expected to lead to a much larger radius 
parameter along the direction of the total momentum of detected two particles than 
that perpendicular to both this direction and the beam direction. 

From the emission function, which is given by the space-time and momentum 
distribution of particles at freezeout in the AMPT model, we can evaluate the 
correlation function C(Q,K) of two identical pions including their Coulomb inter- 
action. Here, Q and K are the relative and total momenta of the two pions. The 
six-dimension correlation function is usually shown as a function of the invariant 
momentum (Qi nv = \/—Q 2 ) or as a function of the projection of the relative mo- 
mentum Q in the "out-side-long" (osl) system [||[ |3^| , defined by the beam direction 
(Qiong), the direction along the total transverse momentum of the two pions (Qout), 
and the direction orthogonal to the above two directions (Qsidc)- 

In the left panel of Fig. ||, we show the one-dimensional projections of the 
correlation function for midrapidity pions with transverse momentum 125 < p±_ < 
225 MeV/c in Au+Au collisions at ^fs = 130A GeV. Also shown is the measured 
7T _ correlation function by the STAR collaboration without correcting the effect 
due to Coulomb interaction [ ^3. In evaluating the one-dimensional projection of 
the correlation function onto one of the Qout, Qsidc, Qiong axes, we have integrated 
the other two Q components over the range — 35 MeV/c. Dash-dotted curves 
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Fig. 4. Left panelrCorrelation functions for midrapidity pions from central Au+Au 
collisions at s/s = 130 A GeV. Coulomb- uncorrected correlation functions from the 
STAR collaboration for low pt tt~ at midrapidity from central collisions are shown 
by filled circles. Dash-dotted curves are results from the default AMPT model, while 
other curves are from the extended AMPT model with string melting and various 
values for the parton scattering cross section. Right panel: Distribution of the 
emission source from the AMPT model for midrapidity pions with 125 < p±_ < 225 
MeV/c in the x out — £ s ide (upper panel) and a; out — t (lower panel) spaces. 



in the figure are from the default AMPT model with a parton scattering cross 
section of a p = 3 mb, while other curves are from the extended AMPT model with 
string melting but different values of parton scattering cross section. It is seen 
that with string melting both the width of the Q lnv correlation function and its 
height decrease with increasing er p . To reproduce the measured one-dimensional 
correlation functions from the STAR collaboration, we need a parton scattering 
cross section of about 10 mb in the extended AMPT model with string melting. 

In the right panel of Fig. ^, we show the distribution of the emission source 
from the AMPT model for midrapidity pions with 125 < p± < 225 MeV/c in the 
x ut — £sidc and Xout — t spaces. It is seen that the emission source shows a large 
halo around a central core. The halo consists of not only pions from decays of 
long-lived resonances such as the u> but also thermal pions. The latter contribution 
to the halo becomes increasingly important when the collective expansion velocity 
of the source is large as a result of increasing parton cross section. The emission 
source also shows a strong correlation between x ou t and t, with the width of its x ou t 
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distribution increasing with the emission time t. 
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Fig. 5. Left panel: Radii i?o Ut , ^sidc anc ^ -^long an< ^ emission duration At of 
the emission source for midrapidity pions with 125 < pj_ < 225 MeV/c in central 
Au+Au collisions at s/s — 130 A GeV as functions of the average transverse mass 
of the pion pair from the AMPT model with string melting and parton scattering 
cross section a p = 15 mb. Right panel: Same as left panel for the radius parame- 
ters obtained from fitting the two-pion correlation by a Gaussion function in their 
relative momentum. 

The size of emission source can be estimated from the curvature of the correla- 
tion function at Q = 0, i.e, the second moments of the emission function from the 
AMPT model. In the left panel of Fig. [|, we show the dependence of the source 
radii of midrapidity pions (—0.5 < y < 0.5) on the average transverse mass tot of 
the pion pair from the AMPT model with string melting and a parton cross section 
of 10 mb. At low transverse mass, all three radii i?o Ut , -R^dei an( i R\ong are l ar g e 
and have values between 10 and 25 fm/c. The emission duration At also lasts a 
long time of more than 25 fm/c. As the average pion transverse mass increases, 
both the emission duration and the source radii, particularly R\ ongl are found to 
decrease. The ratio -Ro U t/^side ^ as a vauie °^ a bout 1.2 for most values of tot and 
is smaller than the prediction based on the hydrodynamical model with freezeout 
treated via the hadronic transport model [ [35| . 

Empirically, radii of the emission source are obtained from fitting the measured 
correlation function by a Gaussian function in the relative momentum of the pion 
pair. Results from such a fit to the correlation function obtained from the AMPT 
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model with string melting and a p = 10 mb are shown in the right panel of Fig. 
^ as functions of the average pion transverse mass. The fitted radius parameters 
Rout, fisidc, arid i?iong are about a factor of two to three smaller than correspond- 
ing source radii evaluated directly from the emission function, indicating that the 
emission source from the AMPT model deviates appreciably from a Gaussian one. 
Furthermore, the ratio -R ut/-Rsidc of the fitted Gaussian radius parameters is much 
closer to one than that from the source radii obtained from the emission function. 

6. Conclusions 

Using a multiphase transport model, which includes both initial partonic and final 
hadronic interactions, we have studied the rapidity distributions of charged particles 
such as protons, antiprotons, pions, and kaons in heavy ion collisions at RHIC. 
With the model parameters constrained by data from central Pb+Pb collisions at 
^/s = 17 A GeV at SPS, theoretical results on the total charged particle multiplicity 
at midrapidity in central Au+Au collisions at ^fs = 56^4 and 13(M GeV agree quite 
well with the data from the PHOBOS collaboration. Also, the predicted p/p ratio 
is consistent with that from the STAR and BRAHMS collaborations. Furthermore, 
the centrality dependence of the charged particle rapidity distribution measured by 
the BRAHMS collaboration is well reproduced. These hadronic obscrvables are, 
however, less sensitive to the initial partonic interactions than to the final hadronic 
interactions. On the other hand, partonic effects are important for both the elliptic 
flow of charged particles and the two-pion correlation function. To reproduce the 
measured elliptic flow and two-pion correlation function by the STAR collaboration, 
we need not only to convert the initial strings to partons but also to use a large 
parton-parton scattering cross section. Our results thus demonstrate the possibility 
of studying the dynamics of partonic matter in heavy ion collisions at RHIC. To 
further test these partonic effects, it is necessary to study their dependence on the 
initial parton distribution, the inelastic scatterings in the partonic stage, and the 
mechanism for melting the initial strings as well as for hadronization. Such an 
improved study based on the AMPT model will be very useful in finding the signals 
for the quark-gluon plasma formed in heavy ion collisions at RHIC. 
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